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ABSTRACT 

Cosmic ray protons interacting with gas at the mean density of the interstellar medium in starburst galax- 
ies lose energy rapidly via inelastic collisions with ambient nuclei. The resulting pions produce secondary 
electrons and positrons, high-energy neutrinos, and 7-ray photons. We estimate the cumulative 7-ray emission 
from starburst galaxies. We find a total integrated background above 100 MeV of w 10"^ GeV cm"^ s"' sr"' 
and a corresponding specific intensity at GeV energies of vly « 10"^ GeV cm"^ s~' sr"^ Starbursts may thus 
account for a significant fraction of the extra-galactic 7-ray background. We show that the FIR-radio correla- 
tion provides a strong constraint on the 7-ray emission from starburst galaxies because pions decay into both 
7-rays and radio-emitting electron/positron pairs. We identify several nearby systems where the potential for 
observing 7-ray emission is the most favorable (M82, NGC 253, & IC 342), predict their fluxes, and predict a 
linear FIR-7-ray correlation for the densest starbursts. If established, the FIR-7-ray correlation would provide 
strong evidence for the "calorimeter" theory of the FIR-radio correlation and would imply that cosmic rays in 
starburst galaxies interact with gas at approximately the mean density of the interstellar medium (ISM), thereby 
providing an important constraint on the physics of the ISM in starbursts. 

Subject headings: galaxies:starburst — gamma rays:theory, observations — cosmology:diffuse radiation — 
ISMxosmic rays — radiation mechanisms:non-thermal 



1. INTRODUCTION 

The magnitude of the extra-galactic 7-ray background is 
uncertain, primarily due to the presence of foreground con- 
taminants (compare Keshet et al. 2004, Strong et al. 2004, 
Sreekumar et al. 1998). GLAST^ will have an order of 
magnitude better sensitivity at GeV energies than previous 
experiments and should provide important constraints on 
the unresolved extra-galactic 7-ray emission. A number of 
potential contributors to the GeV 7-ray background have 
been discussed in the literature including blazars (Stecker & 
Salamon 1996), galaxy clusters and groups (Colafrancesco 
& Blasi 1998; Dar & Shaviv 1995), intergalactic shocks 
and structure formation (Loeb & Waxman 2000; Keshet et 
! al. 2003; Miniati 2002), dark matter annihilation (e.g., UUio 
et al. 2002; Elsasser & Mannheim 2005), and normal star- 
forming galaxies (Pavlidou & Fields 2002). In this paper we 
calculate the contribution of starburst galaxies to the GeV 7- 
ray background. In particular, we assess the possibility that 
starbursts are cosmic ray proton calorimeters: the dense ISM 
of these systems acts as a beam dump for the total energy in- 
jected by supernovae into cosmic ray protons, with a portion 
of the proton energy emerging as 7-rays and high-energy neu- 
trinos. 

In ^we review Volk's (1989) electron calorimeter model 
for the observed FIR-radio correlation (see Thompson et 
al. 2006, hereafter T06). We then argue that starbursts may 
also be proton calorimeters. In |3l we predict the 7-ray flux 
from starburst galaxies and highlight which nearby systems 
are most Ukely to have detectable 7-ray fluxes. We argue that 
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the observed FIR-radio correlation provides important con- 
straints on the 7-ray emission from starbursts because pions 
from inelastic proton-proton collisions produce both 7-rays 
and secondary electrons andpositrons, which then produce 
radio synchrotron. Section then discusses the cumulative 
starburst contribution to the extra-galactic 7-ray background. 
Estimates similar to those presented here for the high-energy 
neutrino background from starbursts have recently been made 
by Loeb & Waxman (2006). 

Although the 7-ray flux from individual starbursts has 
been estimated by several authors (e.g., Torres 2004, Tor- 
res et al. 2004, Cilhs et al. 2005, PagHone et al. 1996, and 
Blom et al. 1999), this paper estimates the expected extra- 
galactic background, defines the requisite conditions for pro- 
ton calorimetry, and explicitly connects the 7-ray predictions 
with constraints from the FIR-radio correlation. 

2. STARBURST GALAXIES AS COSMIC RAY CALORIMETERS 

Observations of both normal star-forming and starburst 
galaxies reveal a linear correlation between their radio (pri- 
marily synchrotron) and FIR luminosities (van der Kruit 1 97 1 , 
1973; Helou et al. 1985; Condon 1992; Yun et al. 2001; Mur- 
gia et al. 2005; Murphy et al. 2006). Volk (1989) argued 
that the FIR-radio correlation arises because the synchrotron 
cooling time for relativistic electrons generated in supernova 
shocks is significantly shorter than their escape time from the 
galactic disk (the "calorimeter theory"). In this limit, the cos- 
mic ray electrons radiate nearly all of the energy supplied to 
them by supernovae. The calorimeter theory is a particularly 
compelling explanation for the FIR-radio correlation in star- 
burst galaxies, where the synchrotron and inverse Compton 
(IC) cooling times of cosmic ray electrons are much shorter 
than in normal star-forming galaxies (T06). 

A priori, one of the strongest objections to the calorimeter 
theory is that the radio spectra of galaxies are not compatible 
with the expected steepening due to strong synchrotron cool- 
ing (e.g., Condon 1992; Niklas & Beck 1997): in the pres- 
ence of strong cooling, models predict /v « i/"" with awl. 
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while observations find a ~ 0.75, even in starbursts (e.g., 
Condon et al. 1991; Condon 1992; Niklas et al. 1997). T06 
argue that if cosmic ray electrons interact with gas at approxi- 
mately the mean density of the ISM, ionization and relativistic 
bremsstrahlung losses will systematically flatten the nonther- 
mal spectra of starburst galaxies and simultaneously maintain 
the linearity of the FIR-radio correlation. This strongly sup- 
ports Volk's calorimeter model for the FIR-radio correlation. 
In the model of T06, the conclusion that cosmic rays interact 
with gas at roughly the mean density of the ISM is required to 
account for the radio spectra of starbursts in the presence of 
strong synchrotron cooling. This requirement has an impor- 
tant consequence for cosmic ray protons, accelerated together 
with electrons in supernova remnants. 

If cosmic ray protons interact with gas at the mean den- 
sity of the ISM in starbursts, then they lose energy rapidly via 
inelastic proton-proton collisions with ambient nuclei. This 
interaction generates both charged and neutral pions, which 
subsequently produce secondary electrons and positrons, high 
energy neutrinos, and photons, respectively. For protons of 
GeV energy, the total inelastic proton-proton cross section is 
Ri30mb (with inelasticity '^0.5; e.g., Gaisser 1990) so that the 
cosmic ray proton energy-loss timescale is 

T„, « 7 X 10^1-1 yr, (1) 

where « is the number density in units of cm"^. The energy- 
loss timescale decreases by less than a factor of two over six 
orders of magnitude in proton energy. 

If the pion production timescale within a galaxy is shorter 
than the escape timescale of relativistic protons, then the pro- 
tons will lose most of their energy to secondary photons, elec- 
tron/positron pairs, and neutrinos before escaping the galaxy. 
The galaxy would then be a cosmic ray proton calorimeter 

The escape timescale for cosmic ray protons is uncertain. 
Both energy-dependent diffusion and advection in a galactic 
outflow likely contribute. In the Galaxy, the diffusive escape 
timescale for 10 GeV protons is inferred to be Tdift ~ 3 x 10^ 
yr near the solar circle (Garcia-Munoz et al. 1977; Connell 
1998). Such protons interact with a mean total column density 
of '^lO g cm"^ before escaping, and the column decreases at 
higher energies roughly as E"'/^ (e.g., Engelmann et al. 1990; 
Webber at al. 2003). These observations imply that 10 GeV 
cosmic rays interact with matter of density n ^ 0.2 cm""*, 
somewhat below the average density of the ISM in the gas 
disk of the Milky Way. Using equation JO we infer that 10 
GeV protons lose roughly 10% of their energy to pion produc- 
tion before escaping. Assuming that diffusive losses proceed 
similarly in other star-forming galaxies, these results imply 
that there is a critical gas surface density I^crn above which 
protons will lose most of their energy to pion production. Tak- 
ing the Milky Way's surface density to be Smwq ~ 2.5 x lO"-' 
g cm~^ at the solar circle (e.g., Boulares & Cox 1990) we es- 
timate that 

- 10 EMWe ~ 0.03 g cm-l (2) 

For starburst galaxies, the assumption that diffusion dom- 
inates the escape of cosmic rays likely overestimates the es- 
cape timescale. In these systems escape is probably domi- 
nated by advection out of the gas mid-plane in a starburst- 
driven wind (as, e.g., in M82; see Klein et al. 1988, Seaquist 
& Odegard 1991). An estimate for the escape timescale is 
then 

Twind ~ h/V « 3 x lO^/iiooVjoo yr, (3) 



where /z loo = /j/ 1 00 pc is the gas scale height and V300 = V /300 
km s"' is a typical wind velocity (e.g., Martin 1999). The con- 
dition Tpp < Twind is then the criterion for cosmic ray protons 
to lose a significant fraction of their energy to pion produc- 
tion. This criterion may also be written in terms of a critical 
gas surface density; 

I]™'^«0.3y3oogcm-2. (4) 

Although uncertain, equations (|2j and (0} imply that galax- 
ies with > 0.03 - 0.3 g cm~^ are likely to be proton 
calorimeters. A number of local starburst galaxies satisfy this 
criterion (see Table 1). 

3. 7-RAY EMISSION FROM STARBURST GALAXIES & 
THE FIR-7-RAY CORRELATION 

To estimate the 7-ray emission associated with pion pro- 
duction in starbursts, we assume that the star formation rate 
(M*) is related to the total IR luminosity (Ltir [8- 1000] ^m; 
see, e.g., Calzetti et al. 2000, Dale et al. 2001) by Ljir = eM*^^ 
(e is an IMF-dependent constant), that the supernova rate per 
unit star formation, Fsn, is a constant fraction of Mi,, and that 
5?7o.o5% of the energy of each supernova (£51 =£'sn/10^' ergs) 
is supplied to relativistic protons. With these assumptions, the 
total energy per unit time injected in cosmic ray protons is 

LcR„ ~ 4.6 X lO'^iR^si 7?o.05 Pn, (5) 

where (3n = (rsN/e)/17MQ ; for continuous star formation 
over ~10** yrs, this ratio has only a very weak dependence 
on the assumed IMF, since high mass stars dominate both the 
total luminosity and the supernova rate. In addition, for star- 
formation timescales of 3 x 10^-10^ yrs, f3 increases by only 
a factor of approximately 1.5, from «15 to «23 (e.g., Lei- 
therer et al. 1999). We further assume that starburst galaxies 
are proton calorimeters (Tpp < Tgsc, 'Sg > ^cnt) and — for the 
purposes of an analytic estimate — that the 7-ray spectrum is 
flat at and above GeV energies; this requires that cosmic ray 
protons are injected with a spectral index /:> « 2 (e.g.. Bland- 
ford & Eichler 1987). Numerical calculations for different p 
are shown in Figure[nand discussed below. 

Because the cross-sections for neutral- and charged-pion 
production are similar, w 1/3 of the proton energy goes into 
neutral pions, which subsequently decay into 7-rays. Thus, 
under the assumption of proton calorimetry, the expected 
7-ray luminosity can be directly related to Lcr,, and Ljir 
(eq. ||5l). For a flat spectrum and for energies >GeV, we find 
that 

i.L,(GeV) « 1 .5 X 10-4^^^^ « lO'^ 05 l^ir, (6) 
ln(7mflx) 

where ''yma.x{= 10^) is the assumed maximum Lorentz factor of 
the accelerated protons and we have dropped the dependence 
here and below on the ratio (3 = (Fsn / e). 

Approximately 2/3 of the energy in cosmic ray protons 
goes to charged pions, which subsequently decay into high 
energy neutrinos and electron/positron pairs. The pairs re- 
ceive « 1/4 of the pion energy (Schlickeiser 2002). Thus 
« 0.8770.05% of the energy of each supernova goes into sec- 
ondary electrons and positrons. Assuming that the syn- 
chrotron cooling time of the secondary pairs is less than their 
escape time (Volk 1989; T06), we find a FIR-radio correlation 
of the form 

j/L^(GHz) w7.5 X io-s ^?o»5^TiR ^ ^ lO-So.osixiR- (7) 
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The observed correlation is uL^(lAGliz) « 1.1 x 10~''Ltir 
(Yun et al. 2001).^ For rj = 0.05, equation |7) somewhat 
over-predicts the observed FIR-radio correlation in starbursts. 
However, this prediction neglects ionization, inverse Comp- 
ton, and bremsstrahlung losses, which may be competetive 
with synchrotron losses in starbursts and will reduce the mag- 
nitude of vLj^iGHz) / LjirJ For example, in the models of 
T06, which simultaneously explain the linearity of the FIR- 
radio correlation and the radio spectra of starbursts (their 
Figs. 3&4), bremsstrahlung and ionization losses account 
for «50% of the energy lost by cosmic ray electrons. Note 
that even if inverse Compton and bremsstrahlung losses com- 
pletely dominate ionization and synchrotron losses for sec- 
ondary electrons/positrons, their contribution to the total X- 
ray and 7-ray emission from starbursts is at most 1/2 of the 
contribution from neutral-pion decay; inverse Compton and 
bremsstrahlung will, however, dominate the emission below 
-100 MeV. 

Equations (|6j and (0 highlight the importance of the FIR- 
radio correlation for constraining the 7-ray emission from 
starbursts. Statistically, the 7-ray emission from galaxies 
cannot exceed that predicted by equation (|6j without violat- 
ing the radio constraints implied by the FIR-radio correla- 
tion (eq. 0. The only way out of this conclusion is if ion- 
ization, inverse Compton, and bremsstrahlung losses signif- 
icantly exceed synchrotron losses for cosmic ray electrons 
and positrons. Although possible, it would be remarkable 
if synchrotron losses were highly sub-dominant in starbursts 
and yet somehow the linearity of the FIR-radio correlation is 
maintained, both within different starburst galaxies and be- 
tween starbursts and normal star-forming galaxies. 

Equation Q is only consistent with the observed FIR-radio 
correlation for rj — 0.05, i.e., if —5% of the supernova energy 
is supplied to relativistic protons. Although this efficiency 
of cosmic ray production is reasonably consistent with the 
energetics of galactic cosmic rays (e.g., Schlickeiser 2002) 
and with the 7-ray luminosity of the Milky Way (Hunter et 
al. 1987), some models of non-linear diffusive shock acceler- 
ation imply 77 ~ 0.5 (e.g., Ellison & Eichler 1984; Ellison et 
al. 2004).^ If indeed 77 > 0.05, then the FIR-radio correlation 
implies that either starburst galaxies are not proton calorime- 
ters or that synchrotron losses are highly sub-dominant with 
respect to ionization, inverse Compton, and bremsstrahlung 
losses. 

Equations (|6j and Q imply that the 7-ray fluxes of the 
densest starbursts should be linearly proportional to their 
radio and IR fluxes in the proton calorimeter limit. Ta- 
ble n provides a sample of galaxies chosen to have large 
IR fluxes and gas densities (see also Paglione et al. 1996; 
Blom et al. 1999; Torres et al. 2004; Torres 2004; Domingo- 
Santamarfa& Torres 2005; Torres & Domingo-Santamaria 
2005). For the latter, we have applied a cut > 0.02 g cm"-, 
plausibly in the proton calorimeter limit; the implied mean 

' LpiR, as defined in Yun et al. (2001) from the 60^tm and lOO/xm IRAS 
fluxes (Helou et al. 1988), is generally a factor of Ri 1.75 smaller than Ltir 
for starburst galaxies (see Calzetti et al. 2000). 

' Equation Q also neglects the contribution to the radio emission from 
primary electrons accelerated directly in supernova shocks. Typical esti- 
mates for the ratio of the total energy injected into protons and primary elec- 
trons are ~20 — 40 (e.g., Beck & Krause 2005), in which case secondary 
electrons/positrons likely dominate the radio emission in starbursts (see T06; 
Rengarajan 2005). 

* However, Volk et al. (2003) argue that standard spherical models of non- 
linear shock acceleration over-estimate the efficiency of cosmic ray produc- 
tion by a factor of ~5. 



gas densities range from ^^20 cm"-' for IC 342 to —10'* cm"' 
for Arp 220. The systems in Table ^ are listed in order of 
decreasing TIR flux density, which corresponds to decreasing 
7-ray flux in the proton calorimeter limit (eq. (6)). Several 
systems just miss our cut on Sg, including the nucleus of NGC 
5236, NGC 4631, and NGC 3521.'' Because typical errors in 
Eg are of order ±0.2 dex (e.g., Kennicutt 1998) and because 
of uncertainty in J^„it (©, these systems may also prove to 
be bright sources in the 7-ray sky. 

The 7-ray fluxes predicted by equation (|6j at 1 GeV from 
Ltir are also provided in Table For comparison, the ex- 
pected sensitivity of GLAST at 1 GeV for a 1 year integration 
is « 2 X 10""' GeV cm"2 s"' (see Fig.0.'° Table [T] shows 
that there are 3 nearby starbursts above this limit (M82, NGC 
253, & IC 342) and approximately 10 other systems with pre- 
dicted fluxes within a factor of 4 of GLAST's detection limit. 
We note that NGC 3690 and Arp 220 would be particularly 
important for testing our predicted FIR-7-ray correlation over 
a large dynamic range in Ljir. 

Figure^shows the 7-ray spectrum, as calculated using the 
formalism of Aharonian & Atoyan (2000) (their §4.1), as a 
function of energy for the specific cases of NGC 253 {dashed 
lines) and Arp 220 {dotted lines), for p = 2.0-2.4, where p 
is the slope of the proton injection spectrum. We first specify 
p, and then we calculate the spectrum from equation (22) of 
Aharonian & Atoyan (2000), adopting their fit for the cross 
section and their value of /T^ w 0.17 (see also the Appendix 
of Domingo-Santamaria & Torres 2005). We normalize the 
resulting spectrum, such that the integral over all 7-ray ener- 
gies is set by the second equality in equation (|6}. Comparing 
our analytic estimate in equation (|6j with this numerical cal- 
culation, we find that a nearly flat p = 2.1 proton spectrum 
gives the best agreement at GeV energies. 

In order to be consistent with the FIR-radio correlation 
as given in the second equality in equation Q, all calcu- 
lations with p 2.1 are normalized to the p = 2.1 spec- 
trum at 1 GeV. Note that this constraint implies that p and 
77 are not independent. A constant flux at 1 GeV requires 
{r],p) « (0.09,2.0), (0.05,2.1), (0.03,2.2), (0.025,2.3), and 
(0.02,2.4). Because the average magnetic field strengths in 
starbursts may range from Z? — 0.1 mG (e.g., NGC 253 or 
M82) to B - 10 mG (e.g., Arp 220) (see T06 for a discussion), 
and because the energy of GHz-emitting electrons/positrons is 
proportional to Z?"'/^, our normalization of the p =/2.1 curves 
at 1 GeV (rather than the proper normalization of the 1 GHz 
radio emission) introduces a factor of w 2 uncertainty in the 
predicted 7-ray emission. 

As Figure[T]emphasizes, NGC 253 (and M82, see TableQ 
should be readily detected with GLAST. Existing observations 
by HESS at TeV energies (Aharonian et al. 2005) yield an up- 
per limit of < 1 .9 X 10"'^ ph cm"^ s"' sr"' (for a point source; 

> 300 GeV) for NGC 253 that is very close to our predic- 
tion for p = 2, suggesting that longer integrations with HESS 
should be able to detect NGC 253 or significantly constrain 
the high energy spectral slope (see also Domingo-Santamaria 
& Torres 2005; Torres & Domingo-Santamaria 2005). 

In two cases — NGC 253 and Arp 220 — we can compare 

' The Seyfert 2 galaxy NGC 1068 with Eg 0.1 g cm^^ (Heifer et 
al. 2003) would be fourth in TableQwith uF^iGeY) ^ 11.8/ X lO"" GeV 
cm"^ s"' where / is the fraction of the infrared emission due to the nuclear 
starburst rather than the AGN. 

See http://people.roma 2.infn.it/^glast/,^ for a comparison of 7-ray tele- 
scope sensitivities. 
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7-Ray Background (Strong et al. 2004) 
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Fig. 1 . — 7-ray flux and specific intensity as a function of 7-ray energy. Tlie data points with error bars show the extra-galactic 7-ray background as inferred 
by Strong et al. (2004) from the EGRET data. The open square shows the upper limit on the background derived by Keshet et al. (2004) at 1 GeV. The thick solid 
lines are the point source sensitivities of EGRET and GLAST for a one year all-sky survey and the HESS sensitivity curve for a Scr point source detection using a 
50 hour integration. The soHd lines show the total starburst background calculated with p = 2.0 — 2.4 and assumin g ri, ,, = 0.3, Ha = 0.7, Ho=l\ km s"' Mpc"' , 
/(z) = min[l,0.9z4-0.1], and p*(z) from RSF2 in Porciani & Madau (2001). For this model f^.^\ ^ 0.8 (see eg llOt . The expected signal from the starbursts 
NGC 253 (clashed lines) and Arp 220 (dotted lines) with p = 2.0-2.4 (see ^ eqs. |S| & Q) are also shown. For other starburst galaxies potentially visible by 
GLAST and HESS, see TableQ For both the 7-ray background and the individual starbursts, the p = 2.1 curves are normalized such that the total integrated 7-ray 
luminosity is 1.5 X lO^ryo.os times the total integrated stellar luminosity (eq. |S|). Because of the constraints on the radio emission from stai'bursts, all p^2.\ 
curves are normalized to the p = 2.1 calculation at 1 GeV, which implies that in our models, p and 77 are not independent (see discussion in Note that the 
effects of attenuation of the background spectrum above ^ 10 GeV have been neglected (see, e.g., Madau & Phinney 1996). 



our estimate for the 7-ray luminosity with more detailed cal- 
culations by Torres (2004) and Domingo-Santamarfa& Torres 
(2005). To make an apposite comparison, it is necessary to 
scale our results to their assumed paramaters. Using = 0. 165 
and p = 2.2, Torres (2004) finds for Arp 220 that j//v(GeV) « 
3 X 10"'" GeV cm"^ s"'. Rescaling our 7/ « 0.03 and p = 2.2 
result, and adjusting for our somewhat (« 20%) higher TIR 
flux, we find that i//v(GeV) « 2.8 x 10""' GeV cm'^ s"'. 
For NGC 253, Domingo-Santamaria& Torres (2005) esti- 
mate ;//v(GeV) sa 2 X 10"'' GeV cm"^ s"' with = 0.1 and 
p = 2.2. Again scaling our results, we find that j//v(GeV) « 
2.3 X 10"^ GeV cm"^ s"' for their assumed parameters. The 
correspondence between our predicted fluxes and these more 
detailed calculations is thus excellent. These studies also 
find, in agreement with our assumptions, that tt" decay dom- 
inates 7-ray emission at GeV energies, that secondary elec- 
trons and positrons dominate the radio emission at GHz fre- 
quencies, and that both systems are safely in the proton and 



electron/positron calorimeter limits; the cooling timescale is 
less than the escape timescale. 

The predicted fluxes for Arp 220 and NGC 253 above - for 
Torres (2004) and Domingo-Santaman'a & Torres's (2005) pa- 
rameters - are a factor of w 5 and 3 times higher than those 
given in Table (primarily because of the larger values of 
77). Given the constraints imposed by the FlR-radio correla- 
tion, it is therefore perhaps surprising that Torres (2004) and 
Domingo-Santamarfa& Torres's (2005) are able to fit the ra- 
dio emission from Arp 220 and NGC 253 with the signifi- 
cantly larger values of rj they assumed. As per the discus- 
sion after equation Q, this is only possible in the proton 
calorimeter limit if ionization, bremsstrahlung, and/or inverse 
Compton losses dominate synchrotron losses. Indeed, for the 
individual nuclei of Arp 220 and for the inner starburst of 
NGC 253, Torres (2004) and Domingo-Santamaria& Torres's 
(2005) parameters (gas densities, magnetic field strengths, 
etc.) imply that the ratio of the synchrotron cooling timescale 



to the total electron energy loss timescale (including all pro- 
cesses; dominated by ionization losses) is ^5. Thus, their 
results are fully consistent with the arguments of our paper. 
However, we believe that their predicted 7-ray fluxes are op- 
timistic because, as noted above, the tightness and linearity of 
the FIR-radio correlation are difficult to explain if synchrotron 
losses are highly sub-dominant in aU star-forming galaxies. 

4. THE STARBURST 7-RAY BACKGROUND 

The total diffuse extra-galactic TIR background 
from star formation is typically estimated to be 
= FTiR/20nW m-2 sr-'w 1, with roughly half of 
the contribution to the background coming from z > 1 (e.g., 
Nagamine et al. 2006). Equation (|6j for the 7-ray flux from a 
given starburst in the proton calorimeter limit then implies a 
total integrated 7-ray background of 

F^^lAx 10-^0.05/0.75^-2™ GeV S-' cm'^ sr"' (8) 

where /0.75 = /cai/0.75 and 

_ Mz)f(z)dz , f r Mz)dz \ 

is the fraction of the TIR extragalactic background produced 
by starburst galaxies that are in the proton calorimeter limit 
(i.e., Ej > Scrit; f{z) is that fraction at every z, is the 
comoving star formation rate density, and£'(z)= [ri„,(l +z)-' + 
J7a]'''^- For a flat 7-ray spectrum, equation implies a spe- 
cific intensity at GeV energies of 

z//,(GeV) « 10-'7]o.05./o.75^'2o'' GeV s"' cm'^ sr"' . (10) 

The estimated 7-ray background depends sensitively on the 
fraction of star formation that occurs in the proton calorime- 
ter limit, i.e., on f(z) and /cai. We estimate f(z = 0) ^ 0.1 
from the fraction of the local FIR and radio luminosity den- 
sity produced by starbursts (Yun et al. 2001). Confirmation 
that local starbursts like M82 and NGC 253 are 7-ray emitters 
at the level predicted in Table [2 would significantly solidify 
the background predictions of equations (|8} and (I10> by con- 
firming that such starbursts are indeed proton calorimeters. At 
high redshift, a much larger fraction of star formation occurs 
in high surface density systems likely to be calorimeters. For 
example, the strong evolution of the IR luminosity function 
with redshift indicates that is dominated by luminous in- 
frared galaxies at z 1 (Dole et al. 2006); these systems are 
typical of those we expect to be proton calorimeters, in which 
case f(z) ~ 1 at z ^ 1. In addition, the typical gas surface 
density of UV selected galaxies, which account for a substan- 
tial fraction of the star formation at high redshift, is roughly 
^ 0.06Mio.5/?g^p(, g cm"^, where Mio 5 is the inferred av- 
erage gas mass in units of 10'°^ M0 and /?6kpc = ^/6 kpc is 
the half-light radius (Erb et al. 2006ab). This surface den- 
sity is comparable to Yicnt estimated in ^ Thus, if we take 
fiz 1) 1 and /9*(z) as in, e.g., Porciani & Madau (2001), 
we find that /cai 1 (eq. 1101 ). with a significant contribution 
to the 7-ray background coming from z « 1 - 2, so that 

/v,oczy/^(GeV)cx/(z«l-2). (11) 

More specifically, taking 57„, = 0.3 and iljy = 0.7, and assum- 
ing a function /(z) = min(0.9z+ 0.1,1), which smoothly inter- 
polates from a small local starburst fraction (/(z « 0) « 0. 1) to 
an order unity starburst fraction at high redshift, we find that 
/cai « 0.8 for models RSFl, RSF2, and RSF3 from Porciani 
& Madau (2001) and for the "Fossil" model of Nagamine et 
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al. (2006). These models also all have FJq^ « 1. As an al- 
ternative, assuming /(z) = min[0.1(l +z)^,0.8], in which the 
starburst fraction increases roughly in proportion to p*(z), 
/cai ~ 0.6 for the same models." 

Using a different methodology, Pavlidou & Fields (2002) 
calculate the contribution to the 7-ray background from star- 
forming galaxies and find a value for i^I^ at GeV energies that 
is 4-5 times higher than that given in equation (I10> . Pavli- 
dou & Fields (2002) did not account for the constraints im- 
posed by the FIR-radio correlation, nor do they explicitly dis- 
cuss the importance of proton calorimetry and the essential 
competition between escape and pion losses (^. By analogy 
with the Galaxy, Pavlidou & Fields (2002) assume that few 
GeV protons lose roughly 10% of their energy before escap- 
ing in all star-forming galaxies in the local universe (z w 0). 
Because roughly 10% of the local TIR luminosity density is 
generated by starburst galaxies, and because we assume that 
all starbursts are calorimeters, local starburst and normal star- 
forming galaxies contribute roughly equally to the 7-ray back- 
ground at z ~ 0. As we have argued above, however, starbursts 
are the dominant mode of star formation at high z and thus 
1-2 Starbursts dominate the contribution of star-forming 
galaxies to the 7-ray background. 

5. DISCUSSION 

Figure ^ shows the numerically calculated 7-ray back- 
ground as a function of 7-ray energy {solid lines), for p = 
2.0-2.4 (see caption for details). The inferred background 
from EGRET from Strong et al. (2004), as well as the con- 
flicting upper limit at 1 GeV derived by Keshet et al. (2004) 
are also shown. 

Our results indicate that 7-ray production from pion decay 
in starburst galaxies contributes significantly to the observed 
background above 100 MeV The considerable uncertainties 
in the 7-ray background determination, primarily because of 
foreground subtraction (see Keshet et al. 2004), complicate a 
more direct comparison. An important point of this paper is 
that the absolute starburst 7-ray flux — both that from individ- 
ual galaxies and that of the background — cannot exceed that 
predicted in Figure ^ and Table 1 without over-producing the 
radio emission from secondary electrons and positrons pro- 
duced in charged pion decay (see ^ eq. Q). This conclu- 
sion is independent of whether starbursts are in fact proton 
calorimeters (as we have assumed), and can only be circum- 
vented if ionization, inverse Compton, and bremsstrahlung 
losses significantly exceed synchrotron losses for cosmic ray 
electrons and positrons in starbursts; such a determination 
would put strong constraints on the origin of the FIR-radio 
correlation (T06; If detected by GLAST, the low energy 
7-ray emission from systems like NGC 253 and M82 (g2| 

" Stecker (2006) presents upper limits on the 7-ray and neutrino back- 
grounds from starbursts that are significantly lower than our estimates here 
(see ^for our neutrino background estimate). There are two reasons for 
this difference. First, Stecker advocates that 22% of the extragalactic TIR 
background is produced at z > 1.2 (see his Table 1) whereas we find that 
all models for p*(z) that are consistent with observations of star formation at 
high-z and the local stellar mass density imply that roughly half of the TIR 
background is produced at z > 1 (e.g., RSFl, RSF2, and RSF3 from Porciani 
& Madau (2001), the "Fossil" model from Nagamine et al. 2006, the dust- 
corrected star formation history given by Schiminovich et al. 2005, and the 
analytic model given by Hemquist & Springel 2003). Second, we argue that 
a fraction of order unity of all star fo rmat ion at z ~ 1 occurs in starbursts 
likely to be calorimeters (our /(z); eg. II II . whereas Stecker quotes 13% in 
the range 0.2 < z < 1-2 and 60% for z > 1.2 (his Table 1). His step function 
model for /(z) significantly underestimates the contribution of z ~ 1 galaxies 
to the 7-ray background. 
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Fig-E Table[3, may help constrain the contribution of inverse 
Compton and bremsstrahlung to electron/positron losses. 

The primary physical uncertainty in our estimate of the 7- 
ray fluxes from starbursts lies in whether cosmic rays do in 
fact interact with gas at approximately the mean density of 
the ISM. In particular, given that galactic winds efficiently 
remove mass and metals from galaxies (e.g., Heckman et 
al. 1990), it is unclear whether the cosmic rays actually in- 
teract with the bulk of the ISM, which is required for pion 
production to be significant Detection of 7-rays from 

starbursts at the level predicted in this paper would thus pro- 
vide an important constraint on the physics of the ISM in 
starbursts and on the coupling between supernovae and the 
dense ISM, which contains most of the mass. In addition, 
the ionization and bremsstrahlung energy-loss timescales for 
electrons/positrons are similar to the proton-proton pion pro- 
duction timescale. Thus, constraints on pion production in 
starbursts via 7-ray observations directly constrain the im- 
portance of ionization and bremsstrahlung losses for shaping 
the radio emission from starburst galaxies (T06). By con- 
trast, detection of (or upper limits on) the 7-ray emission from 
starbursts well below our predictions would rule out the hy- 
pothesis that starbursts are proton calorimeters, although they 
would, in our opinion, still be electron calorimeters. The al- 
ternative possibility suggested by equation (|6|l, that 77 <C 0.05, 
is unlikely, given the energetics of Galactic cosmic ray pro- 
duction (e.g.. Strong et al. 2004). 

We note that the total IR background from the particular 
model of the star formation history of the universe used to 
produce the 7-ray background in Figure^is F^^^ « 2 x 10"^ 
ergs cm"^ s"' sr"', or w 20 nW m"- sr"', with approximately 
80% (/cai ~ 0.8) coming from starburst galaxies. A complete 
calculation, accounting for the contribution from the old stel- 
lar population (which contributes significantly to the back- 



ground at z w 0), finds a total background in starlight that is 
roughly twice as large (e.g., Nagamine et al. 2006). Equation 
Q and, indeed, the existence of the FIR-radio correlation, 
implies that the expected radio background is i//^ (radio) w 
7.5 X 10-i%.o5/[21n(7„,„,,)] « 2.7 x 10"" ergscm'^s"' sr"'. 
A factor of two reduction in (radio) has been included in 
this estimate as per the discussion after equation Q. 

In addition to producing secondary electron/positron pairs 
and 7-rays, pion production is also a significant source of 
high-energy neutrinos, with w 5% of the proton cosmic ray 
energy going into neutrinos. Two-thirds of this energy goes 
to muon-type neutrinos. Therefore, for starburst galaxies that 
are proton calorimeters, we expect a FIR-neutrino correlation 
of the form i^L^(/i- neutrinos) = (2/3)j^L,y(all neutrinos) w 

1 .5 X lO""* 77o.o5iTiR/ Mlmax) ~ 10"^ VO-OsLjiR, identical to the 
7-ray luminosity (eq. |6l). The corresponding cumulative 
extra-galactic /i-neutrino background at GeV energies is then 
the same as given in equation dlOt . Using a very similar esti- 
mate, Loeb & Waxman (2006) have recently argued that star- 
bursts are likely to be an important contributor to the high- 
energy neutrino background. Because the ratio of the 7-ray 
to neutrino luminosity from pion decay depends only on the 
micro-physics of pion production, the 7-ray fluxes from star- 
bursts predicted in this paper will provide a crucial calibration 
of the expected flux of extra-galactic high energy neutrinos. 
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TABLE 1 

Properties of Local Starburst & Star-Forming Galaxies 



Object Name 


D" 
(Mpc) 


(gcm"^) 


(Jy) 


t d 

(Jy) 


(Jy) 


'^1.4GHz'^ 

(mJy) (10 


!/F„(GeV)S 
GeV cm"- 


Data 
s"') Refs. 


NGC3034 (M82)... 


3.3 


0.69 


1313.5 


1355.4 


8302.4 


7657.0 


73 


4,2,3 


NGC 253 


3.5 


0.47 


997.7 


1857.8 


7755.8 


5594.0 


68'' 


1,2,3 


IC 342 


4.4 


0.02 


255.9 


660.7 


2311.6 


2250.0 


20 


1,2,3 


NGC 2 146 


12.6 


0.14 


144.5 


204.2 


1009.8 


1074.6 


8.9 


1,2,3 


NGC 3690 


42.2 


0.04 


119.7 


118.6 


748.0 


678.3 


6.6 


4,2,3 


NGC 1808 


14.2 


0.09 


104.5 


147.2 


729.4 


529.0 


6.4 


1,2,3 


NGC 1365 


18.6 


0.09 


84.1 


185.4 


704.2 


530.0 


6.2 


7,8,3 


Arp 220 


76.6 


12.0 


107.4 


118.3 


691.9 


326.9 


6.1' 


4,2,3 


NGC 891 


7.4 


0.09 


61.1 


198.6 


623.4 


701.0 


5.5 


6,2,3 


NGC 3627 


10.2 


0.04 


56.2 


144.9 


507.3 


458.0 


4.5 


4,2,3 


NGC 660 


12.0 


0.08 


65.5 


103.8 


477.4 


387.0 


4.2 


4,2,3 


NGC 3628 


11.2 


0.02 


48.4 


121.9 


431.9 


525.0 


3.8 


10,3 


NGC 1097 


18.0 


0.10 


46.7 


116.1 


414.0 


415.0 


3.6 


4,2,3 


NGC 3079 


15.9 


3.7 


50.2 


103.5 


407.8 


849.0 


3.6 


4,2,3 



References. — (1) Condon (1987); (2) Kennicutt (1998); (3) Yun, Reddy, & Condon (2001); (4) Condon et al. (1990); (5) Golla & Wielebinski (1994); (6) 
Condon et al. (1996); (7) Sandqvist et al. (1995); (8) Sandqvist (1999); (10) Paglione et al. (2001); (11) Condon et al. (1998) 

"Distance. In most cases, distances computed with H = 11 km s"' Mpc"' . In some cases, distances ai'e taken from the literature. ''Gas surface density. Typical 
uncertainty in this quantity is ±0.2 dex (e.g., Kennicutt 1998). '^IRAS flux density at 60 fim. ''iRAS flux density at 100 /xm. "^Total IR flux density defined 
by logjQ(STiR) = logjQ(2.58560/jm +'^loO/jm) + logio(l-'5), where the factor of 1.75 relates the total [8- 1000]/im flux density to the [40— 120]^m flux density 
(Calzetti et al. 2000; Helou et al. 1988). For any individual system listed here, this bolometric correction is uncertain at the factor of ~1.5 level (see Bell 2003; 
Dale et al. 2001). '^Radio flux density at 1 .4 GHz. ^Predicted 7-ray flux at 1 GeV, scaled from 5tir according to eq. J^, assuming 77 = 0.05 and 'jmax = 10^. ''See 
discussion of comparison with Domingo-Santaman'a & Torres (2005) in ^ 'See discussion of comparison with Torres (2004) in il2l 
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